Anomalous transport induced by nonhyperbolicity 
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In this letter we study how deterministic features presented by a system can be used to perform 
direct transport in a quasi-symmetric potential and weak dissipative system. We show that the 
presence of nonhyperbolic regions around acceleration areas of the phase space plays an important 
role in the acceleration of particles giving rise to direct transport in the system. Such effect can be 
observed for a large interval of the weak asymmetric potential parameter allowing the possibility to 
obtain useful work from unbiased nonequilibrium fluctuation in real systems even in a presence of 
a quasi-symmetric potential. 



Anomalous transport is an emerging field in physics 
and, generally speaking, refers to nonequilibrium pro- 
cesses that cannot be described by using standard meth- 
ods of statistical physics. The investigation of anoma- 
lous transport processes requires a combination of con- 
cept and methods of diverse disciplines, like stochastic 
theory, dynamical systems theory and disordered sys- 
tems [l| . Anomalous transport occurs in a wide realm of 
physical systems ranging from a microscopic level (such 
as conducting electrons) to a macroscopic scale (as in 
global atmospheric events). One of the well- known phe- 
nomena in this category is anomalous diffusion, for which 
the mean-squared-displacement increases with time as a 
power-law i M , where /j^l j|. 

There is a growing interest in anomalous transport 
properties of nonlinear systems presenting nonequilib- 
rium fluctuations, the ratchet systems [3|— 16J . For these 
systems it is possible to surmount the second principle of 
thermodynamics provided we do not have any space or 
time symmetries forbidding direct transport [7J . Ratchet 
systems occur in a variety of physical problems, like uni- 
directional transport in molecular motors 3, Sj, micro 
particles segregation in colloidal solutions]!] L and trans- 
port in quantum and nanoscale systems El ■ 

Previous works have studied ratchet systems with a 
mixed phase space and weak dissipation. For such sys- 
tems it has been shown that the ratchet effect results 
from a connection between weak dissipation, chaotic dif- 
fusion, and ballistic transport due to presence of periodic 
islands 0, HI . The latter, although occurring in the con- 
servative case only, still affects the dynamics in the weak 
dissipative case since island centers become attracting 
fixed points in the weak dissipative case . It has been 
shown that in the absent of islands anomalous transport 
cannot be achieved and the transport currents due to the 
ratchet effect in weakly dissipative systems are related to 
the existence of isoperiodic stable structures 0, [I| . 

In these previous investigations it has been focused the 
strongly asymmetric case, for which the magnitude of po- 



tential asymmetry is comparable with the original poten- 
tial energy of the system. In this letter we show that it is 
not necessary to have such strong asymmetry, in the sense 
that sizeable ratchet currents can be obtained in weakly 
dissipative systems with slightly asymmetric potentials. 
In fact, we claim that the presence of ratchet currents is 
influenced not so much by the potential asymmetry, but 
rather by the existence of strongly nonhyperbolic regions 
in the phase space of weakly dissipative systems. 

As a representative illustration of this effect we con- 
sider a periodically kicked rotor subjected to a harmonic 
potential function, whose dynamics is two-dimensional. 
By a hyperbolic region S we mean a set for which the 
tangent phase space in each point splits continuously 
into stable and an unstable manifolds which are invariant 
under the system dynamics: infinitesimal displacements 
in the stable (unstable) direction decay exponentially as 
time increases forward (backward) [12| . In addition, it 
is required that the angles between the stable and un- 
stable directions are uniformly bounded away from zero. 
Chaotic orbits of dissipative two-dimensional mappings, 
for example, are often nonhyperbolic since the stable and 
unstable manifolds are tangent in infinitely many points. 

The dynamics of a periodically kicked rotor with small 
dissipation and potential asymmetry can be described 
in a cylindrical phase space (— oo x oo) x [0,2-71"), whose 
discrete-time variables p n and x n are respectively the mo- 
mentum and the angular position of the rotor just after 
the nth kick, with the dynamics given by the following 
dissipative asymmetric kicked rotor map (DAKRM) [3] : 

p n+ i = (1 - 7)p„ - K[sm(x n ) + a sin(2x„ + 7r/2)],(l) 
x n +i = x n +p n +i, (2) 

where K is related to the kick strength, < 7 < 1 is 
a dissipation coefficient, and a is the symmetry-breaking 
parameter of the system. The conservative (7 = 0) and 
symmetric (a = 0) limits yield the well-known Chirikov- 
Taylor map In the following we will keep the dissi- 
pation small enough (namely 7 = 2 x 1CP 4 ) in order to 
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FIG. 1. (color online) 2-D histograms (top) and momentum 
probability distributions (bottom) for the DAKRM with 7 = 
0.0002, a = 0.005, and (a,e) K = 6.40; (b,f) 6.92; (c,g) 6.96; 
and (d,h) 7.00. 



highlight the effect of the periodic islands of the conser- 
vative case. Moreover, the asymmetry parameter a will 
be kept small so as to emphasize the role of the nonhy- 
perbolic phase space regions on the anomalous transport. 

The conservative and asymmetric (a ^ 0) case has 
been studied in detail: it has two fixed points (we call 
Pi) given by x R,L — sin -1 ©(a, K) and p R,L — 0, where 

6 = (l - y/l + 8a 2 ± 16™/ K) /4a, which are stable 

centers in the following parameter intervals: < a < 1/4, 
and 6.40 < K < 7.20. These two Pi points are the cen- 
ters of two resonant islands that are actually accelerator 
modes. There are also two (left and right) period-3 fixed 
points (P3) related to secondary resonances around the 
Pi islands. 

In the weak dissipative case, the two Pi points become 
stable foci, their basins of attraction being roughly the 
region occupied by the respective islands (in the conser- 
vative case). Moreover, the chaotic region in the con- 
servative map becomes a chaotic transient in the weakly 
dissipative situation. At K ss 6.92 the right P 3 points 
collides with the right fixed point (p R ,x R ) by a bifurca- 
tion. At the bifurcation point the small attraction basin 
of Pi engulfs the stable manifold of the P3 and turns to 
be accessible to points in a large phase space region. 

The vicinity of the fixed points plays a key role in 
the anomalous transport mechanism, in the same way 
as the islands do for the conservative case. More pre- 
cisely, the wide accessibility of this vicinity near the bi- 
furcation is responsible for large ratchet currents, just as 
the role of the accelerator modes in the non-dissipative 
map. Figs. [UJa-d) depict 2-D histograms for 5000 or- 
bits (each orbit containing 10 3 points) of the DAKRM 
from initial conditions chosen in the phase plane region 
< x < 2tt, — tt < p < n, as well as the corresponding 



momentum probability distributions a(p) (Figs. Qle-h)). 
For K = 6.40 there is a quasi-symmetric situation, the 
neighborhood of the two Pi fixed points (left and right) 
being seldom visited [Fig. (Ha)]. Since the left (right) re- 
gion is responsible for a negative (positive) increase of the 
momentum transport, we observe that for this parame- 
ter value the momentum distribution function is nearly 
symmetric, with a Gaussian shape [Fig. Q3e)], resulting 
in a null transport. 

Symmetry-breaking effects start to be noticeable after 
K = 6.40 and turn to be maximum in the bifurcation 
at K w 6.92, where only the vicinity of the left Pi is 
scarcely visited by orbits of the map [Fig. Hfb)]. This 
effect is triggered by the bifurcation whereby the right 
P3 fixed point collides with the right Pi point and turns 
its vicinity easily accessible (not shielded), what is re- 
flected in the asymmetric right tail in the momentum 
distribution function [Fig. QJf)]. The vicinity of the left 
Pi is not yet affected since the collision process did not 
occur yet for left Pi and P3. As a result, a net trans- 
port current is generated. However, this is more an ef- 
fect of the bifurcation (due to nonhyperbolicity) than of 
the symmetry-breaking itself. In other words, if there 
is weak symmetry-breaking but no bifurcation (and no 
shield process), the ratchet effect will not occur, at least 
with the magnitude we observed in this example. 

Not too far from the bifurcation, (K = 6.96), the 
vicinities of both Pi fixed points become now regularly 
visited [Fig. W[c)] since positive and negative currents, 
caused by the right and left vicinity is counterbalanced. 
For this case the momentum distribution function is 
again approximately symmetric [Fig. QJg)] but present- 
ing right and left tails. The situation changes again af- 
ter the left Pi and P3 fixed points collide [Fig. []Jd)], 
through the same bifurcation mechanism described for 
their right counterparts. The increase in the accessibility 
of the vicinity of the left Pi point leads to an asymmetric 
momentum distribution function [Fig. Qlh)], restoring a 
net transport current. In this last case is noticeable that 
the seldom visited region is very small, nevertheless it is 
enough to inhibit positive currents. Once again the non- 
hyperbolicity of right region seems to be more important 
than the asymmetric situation itself. 

The variation of the average net transport current p 
with the nonlinearity parameter K is depicted in Fig. 
[2ja) for different values of the asymmetry parameter. 
The net current is ensemble-averaged over a large num- 
ber (10 6 ) of initial conditions, each of them being followed 
by a short time (t = 10 3 ) to prevent the system to settle 
down into any fixed point. On varying K we obtain a 
series of positive and negative net transport currents re- 
sulting from the ratchet effect. Curiously the net current 
fluctuates less for both very small and large asymmetry, 
being more sensitive to K for intermediate values of a. 
As we have seen, the appearance of net currents is due to 
the fact that the left and right P3 fixed points (that bifur- 
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FIG. 2. (color online) (a) Ensemble averaged net current for 
different values of the nonlinearity and asymmetry parameters 
of the DAKRM with 7 = 0.0002. (b) Net transport current (in 
colorscale) as a function of the asymmetry and nonlinearity 
parameters. 



cate in pairs for the symmetric case) start to bifurcate at 
different values of K. For the interval 6.40 < K < 7.00, 
corresponding to Fig. Q] the right P3 fixed point bifur- 
cates before the left one, generating a sequence of positive 
and negative net currents with very large peak values. 
For instance, the maximum amplitude of the net trans- 
port current for a = 0.005 is at least three times larger 
than for similar parameters in the case of high asymme- 
try (a = 0.5) [4j. 

By the way of contrast, with a higher asymmetry value 
as those in Ref. |4J such large net currents are observed 
only for larger nonlinearities (8.5 < K < 10), hence they 
are not primarily related to the bifurcations we present. 
In such case the role of the potential asymmetry over- 
comes the bifurcation mechanism presented here. Nev- 
ertheless appreciable net currents can also be acquired 
for asymmetry parameter as small as a — 0.0005 or even 
smaller, although the peak values decrease considerably. 
Such decrease in the net transport current for small val- 
ues of a is expected since for the case of a = no trans- 
port can be observed due to the symmetry of the stan- 
dard map. 

The sensitive dependence of the net transport cur- 
rent on the nonlinearity parameter in the weak asym- 
metry case reminds us of a similar behavior for the 
diffusion coefficient of the conservative and symmetric 
(Chirikov- Taylor) map, caused by the existence of accel- 
erator modes. A more complete comparison is shown in 
Fig. Hlb). If the potential asymmetry is large (a < 0.5) 
sizeable net currents are obtained only for higher values 
of K. This is due to the fact that large asymmetries shift 
the small islands present in the symmetric system phase 
space to some other intervals of K. On the other hand, 
if the asymmetry is very weak (a > 0) we obtain large 
net currents (both positive and negative) for if-values 
considerably lower. The negative transport observed in 
Ref. [H corresponds to the pale cyan area in Fig. E^b) 
for a = 0.5 and around K = 6.5, hence it is not directed 
correlated to the bifurcation mechanism we describe here. 



Moreover, the net current amplitude does not change ap- 
preciably along the blue and red lines in Fig. [Hb) when 
a is varied, emphasizing the role of the nonhyperbolic re- 
gions as an essential requirement for the ratchet effect to 
arise. 

The existence of nonhyperbolic regions in phase space, 
however small they may be, constitutes a deterministic 
mechanism underlying anomalous transport in the pro- 
duction of net currents through a ratchet effect. In or- 
der to quantify the degree of nonhyperbolicity related 
to the phenomena we describe in this letter, let us 
consider an initial condition (po,Xq) and a unit vec- 
tor v, whose temporal evolution is given by v n+ \ = 
J{p n ,x n )v n /\J(p ni x n )v n \, where J{p n ,x n ) is the Ja- 
cobian matrix of the DAKRM. For n large enough, v is 
parallel to the Lyapunov vector u(p, x) associated to the 
maximum Lyapunov exponent \ u of the map orbit be- 
ginning with (j)q,xq). Similarly a backward iteration of 
the same orbit gives us a new vector v n that is parallel 
to the direction s(p,x), the Lyapunov vector associated 
to the minimum Lyapunov exponent A s [Til . HH ] . For re- 
gions where X s < < X u the vectors u(p,x) and s(p, x) 
are tangent to the unstable and stable manifolds, respec- 
tively, of a point (p,x). 

The nonhyperbolic degree of a region S can be studied 
by computing the local angles between the two mani- 
folds 9(p,x) = cos _1 (|« ■ s\), for (p, x) e S [16]. So, 
0(p,x) = denotes a tangency between unstable and 
stable manifolds at (p, x). Let S R ' L — {(p, x) e ft : 
\(p — x) — (p RL ,x R ' L \ < e} be a e-radius neighborhood 
of the right and left Pi fixed points. Results for 9(p 1 x) 
and its distribution function p{9) calculated in both re- 
gions are shown in Fig. [3] for four values of the nonlinear 
parameter K. The dark regions in Fig. [3] correspond to 
strongly nonhyperbolic region surrounding the accelera- 
tion region. 

For K = 6.40, near the fixed points (p R ' L , x R,L ) there 
is a strong nonhyperbolic region which shields the accel- 
eration area [Figs. [H[a-b)]. In fact, almost all 9 values for 
the right area (red curve) and left (blue curve) are con- 
fined in the interval 9 < 7r/16, configuring a strongly non- 
hyperbolic region around both fixed points [Fig. EJc)]. 
By way of contrast, when K = 6.92, only the left fixed 
point neighborhood is shielded, resulting in a large pos- 
itive transport since only the right acceleration region 
is regularly visited [Figs. EJd-e)]. The right Pi and P3 
fixed points suffer a bifurcation and all tangencies of sta- 
ble and unstable manifolds disappear from the right area, 
allowing the trajectories to visit the acceleration region. 
Accordingly Fig. presents different distributions of 
9 values for left and right regions. The red curve (right) 
presents a distribution peak around 9 — tt/8 considerably 
greater than the blue one, leading to absence of shielding 
in the right region. 

In the case of K = 6.96 neither of the areas surround- 
ing the fixed points are shielded by the tangencies of man- 
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for weakly dissipative systems. Let the inner region sur- 
round the P3 fixed point. When the stable and unstable 
manifolds of P3 are almost tangent, the inner region sur- 
rounding P 3 mapped to an outer region tends to zero 
since the angle between the manifolds goes to zero. This 
clearly inhibits transport from inner to outer regions. As 
a result the acceleration region is practically not visited 
and no net transport currents are observed when there 
are nearly parallel manifolds. 

In conclusion we have shown that anomalous trans- 
port displayed by a gitasi-symmetric potential and weakly 
dissipative system is strongly related to the topology of 
the acceleration regions around fixed points displayed 
by the system. The presence of nonhyperbolic regions 
caused by almost parallel unstable and stable manifolds 
can inhibit a chaotic trajectory to visit the neighbor- 
hood of the acceleration region surrounding fixed points 
of the system. This mechanism is closely related to the 
scenario described by the Poincare-Birkhoff theorem in 
area-preserving two-dimensional maps. This dynamical 
phenomenon yields large net transport current in some 
direction even though the potential has an extremely 
small degree of symmetry-breaking. Hence such net cur- 
rents can yield useful work from unbiased nonequilibrium 
fluctuation even with Quasi-symmetric potentials, which 
enlarges the realm of dynamical systems displaying the 
ratchet effect. 

This work is partially supported by CNPq, CAPES 
and Fundagao Araucaria. 



FIG. 3. (color online) 8(p, x) values evaluated from 10 5 initial 
conditions uniformly distributed around (0.2 radius) S^' L and 
the distribution function of 8. 



ifolds, resulting in large positive and negative transport 
currents, but no net current at all [Figs. EJg-i)]. This 
situation, however, is different from the one displayed by 
Figs. [SJa-c), where the regions surrounding both fixed 
points were scarcely visited, hence there is no net current 
since the positive and negative currents are very small. 
Finally, for a higher value of K, only the right accelera- 
tion region is shielded, resulting in a negative transport 
current [Figs. GSj-k)]. The distribution p(9) presents a 
peak near zero for the red curve, confirming the existence 
of a shielded right area [Fig. [3jl)]. 

The source of nonhyperbolicity in these regions is the 
tangencies between stable and unstable manifolds of sad- 
dle orbits embedded in the chaotic region therein. These 
tangencies prevent the system to visit the acceleration 
region. The scenario can be regarded as a counterpart of 
the Poincare-Birkhoff' theorem (that describes the torus 
breakdown of a conservative two-degree of freedom map) 
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